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^ ■ Abstract 

^ \ We consider a simple extension of the Standard Model (SM) that incorporates a Majorana 

pg \ fermion dark matter and a charged scalar particle with a couphng to the SM leptons through 
_il ■ renormalizable terms. Another renormalizable term involving the charged scalar and the Higgs 
O ■ boson gives rise to interactions between the dark matter and SM quarks at the one-loop level, 
. which induce the elastic scatterings between dark matter and nucleus. The same term also 
l; ! affects the effective coupling of the Higgs boson to diphoton through a one-loop diagram with 
the charged scalar. Therefore, our model predicts a correlation between the spin-independent 
cross section for dark matter-nucleus elastic scatterings and a new contribution to the effective 
c^ ■ Higgs boson coupling to diphoton. When the spin- dependent cross section is large enough to 
be tested in future direct dark matter detection experiments, the Higgs-diphoton decay rate 
shows a sizable deviation from the SM prediction. We also consider the case where the fermion 
dark matter is a Dirac particle. Most of discussions is similar to the Majorana case, but we 
find that the magnetic dipole moment of the Dirac fermion dark matter is loop-induced and 
this interaction dominates the spin-independent cross section for dark matter-nucleus elastic 
scatterings. We find that the resultant cross section is about an order of magnitude below the 
current experimental bound and hence can be tested in the near future. 



1 Introduction 

The Wilkinson Microwave Anisotropy Probe (WMAP) observations [1] have estabhshed the 
existence of cold, non-baryonic dark matter (DM), whose abundance is measured to be 

Oj3^^/i2 = 0.112 ±0.006 . (1) 

Since a DM candidate is missing in the standard model (SM), its existence motivates us to 
consider physics beyond the SM. One notable fact about this measured DM abundance is its 
possible relation to physics at the TeV scale. Let us assume that the DM particle is a weakly- 
interacting massive particle (WIMP) that has been in thermal equilibrium with SM particles 
and decouples from the thermal bath when the temperature drops below its mass scale, i.e., 
the DM particle is assumed to be a cold relic. The measured DM abundance, Eq. ([T]), then 
implies that the thermally averaged cross section times velocity for DM annihilations to SM 
particles at the time of decoupling satisfies 

(^^)w. ~ 7^^ ' (2) 

with a coupling constant (g) between the DM and SM particles. Therefore, if the coupling 
constant is of the order of the electroweak gauge couplings, dark matter physics is controlled 
by the TeV scale. 

Motivated by the fact above, we study a simple extension of the SM at the TeV scale that 
incorporates a DM particle and its interaction with SM particles. In this paper, we try to extend 
the SM as minimally as possible. Models for DM have often been associated with solutions 
to the gauge hierarchy problem, a well-known example of which is the neutralino DM scenario 
in minimal supersymmetric standard model. However, the ATLAS and the CMS experiments 
searching for supersymmetric particles at the Large Hadron Collider (LHC) have been reporting 
negative results in all channels, that might cast doubts on the notion of naturalness and the 
validity of the gauge hierarchy problem. In this paper, we adopt the "minimality" as the guiding 
principle for model-building, instead of the naturalness of the electroweak scale. 

The simplest renormalizable extension of SM that incorporates a DM particle is the so-called 
Higgs portal model [2], where a SM gauge singlet scalar is introduced as a DM particle and it 
interacts only with the Higgs boson through a renormalizable term. Unfortunately, this model is 
in tension with XENON 100 direct DM detection experiment [3], except for a limited parameter 
region where the DM mass is nearly a half of the Higgs boson mass and the annihilation of DM 
particles is enhanced by the resonance of s-channel Higgs boson [1] . The reason for the tension 
is that the DM particle and SM quarks interact by the Higgs boson exchange at the tree-level, 
giving rise to too large spin-independent cross section for DM-nucleus elastic scattering. The 
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model evades the bound from the direct DM detection experiments when the DM mass is close 
to a half of the Higgs boson mass and the DM annihilation cross section can be sufficiently 
large for a very small coupling between the DM particle and the Higgs boson |1]. Then it must 
be larger than the Higgs boson mass because otherwise the Higgs boson could decay into DM 
particles with a large branching ratio, which would contradict with the results of the Higgs 
boson searches [5]. When the DM mass is larger than the Higgs boson mass, it is hard to 
observe the signature of DM production at the collider experiments, because DM particles are 
produced only through an off-shell Higgs boson [1]. 

Here we consider the second simplest renormalizable extension of the SM with a fermionic 
DM particle. In this model, we introduce a SM gauge singlet fermion (x) and a scalar (5) 
charged under SM gauge groups. Odd Z2 parity is assigned to S and X; whereas even parity 
is assigned to all the SM particles. Hence, assuming x is lighter than S, the lightest Z2 odd 
particle, x? is stable and a DM candidate. S and x can interact with SM particles through the 
following renormalizable term: 

C D -yi tpiSx + h.c. , (3) 

where ipi denotes a SM particle with flavor index i that has the same SM charge as S. If the 
coupling constant yi is of the order of the electroweak gauge couplings, then the measured DM 
abundance in Eq. ([1]) implies that the particles, x ^i-nd 5", have the mass of order 100 GeV— 1 
TeV. There is a variety of choices for the charge assignment of S to allow the interaction in 
Eq. 03]). In this paper, we consider the case that S is color-singlet and has the mass below 400 
GeV, and the mass difference between 5* and x are smaller than 50 GeV. Such a parameter choice 
is interesting because the charged scalar S has rich phenomenological implications, especially 
on the branching ratio of the Higgs boson into diphoton, but it will take time for the LHC 
experiments to exclude or discover the existence of the scalar and DM particles. 

The heart of this study lies in the fact that since 5* is a scalar particle, no symmetry forbids 
the following renormalizable term: 

C D -k{S^S){H^H) , (4) 

where H is the Higgs boson in the SM. This term induces at the one-loop level an interaction 
between x and SM quarks mediated by the Higgs boson and the spin- independent cross section 
of x-iiucleus elastic scatterings can be accessible by future direct DM detection experiments. 
On the other hand, since S has an electric charge and a mass less than several 100 GeV, the 
same term contributes to the Higgs boson decay to diphoton through a loop diagram, altering 
the Higgs-to-diphoton branching ratio from its SM prediction. Our model thus predicts both a 



signal in direct DM detection experiments and a deviation of the Higgs-to-diphoton branching 
ratio. 

Another important result of this study is that if the DM particle is a Dirac fermion, the 
magnetic dipole moment of the DM particle is induced at the one-loop level. This interaction 
gives a dominant contribution to the DM-nucleus elastic scatterings. We will find that the 
spin-independent cross section for the processes is below the current experimental bounds but 
within the reach of future experiments. 

This paper is organized as follows. In the next section, we define our model. In Sec. 3, we 
investigate DM physics with our model. By calculating the relic abundance of the DM particle 
and comparing our results to the observations, we identify allowed parameter regions of the 
model. We also calculate the spin-independent cross section for the elastic scatterings between 
the DM particle and nucleus mediated by the Higgs boson. In Sec. 4, the effect of the charged 
scalar, S, on the Higgs branching ratio to diphoton is calculated. Then, in Sec. 5, we show 
correlations between the DM-nucleus elastic scattering cross section and the Higgs-to-diphoton 
branching ratio. Finally in Sec. 6, we discuss the case for a Dirac fermion DM particle and 
investigate its phenomenology. 

2 Simple model for Majorana fermion dark matter 

We consider a simple renormalizable extension of the SM which incorporates a fermionic DM 
particle interacting with SM particles. We introduce a SM gauge singlet Majorana fermion (x) 
and a complex scalar field {S) with a hypercharge Y = —1. Odd-parity is assigned to x ^^nd 
S, while even-parity to all the SM fields. The renormalizable terms involving x ^"^^ S are the 
following: 

Centra = '- x^d.x - ^ ru^XX + iD^S)\D,S) - TTilS^ S - Xs {S^ Sf 

- y^ emSx + h.c. - ^{S^S){H^H) , (5) 

where em denotes the SM SU(2)i singlet charged lepton with flavor index i, H denotes the SM 
Higgs doublet field, and yi and k are dimensionless couplings between the new particles and 
the SM fields. The field contents are summarized in Table 1. The DM particle x is assumed to 
be lighter than S, so that x is stable due to the Z2 parity and can be a DM candidate. 

For general values of yj's, the term yiemSx induces flavor-changing neutral current processes, 
which are severely constrained by current experimental results. Here we simply assume that 
yj's are aligned with respect to the SM lepton flavors, i.e., we set yi = y6ei, yS^i or y6ri. 

The masses of x ^-^d S are, respectively, given by m^ and ml + nv"^ /2 = r7i|, after the 
electroweak symmetry breaking by (H) = (f , 0)"^ with v = 174 GeV. Our model is possibly 
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Table 1: The field content of the model 

constrained by the direct slepton search performed by the ATLAS experiment [B]. We here take 
ms > 100 GeV and nis — m^ < 50 GeV, so that the model safely evades the ATLAS bound as 
well as the LEP bound on charged sleptons [7]. 

3 Dark matter physics 

In this section, we constrain the model parameters from the requirement that x be a cold relic 
whose abundance fits with the WMAP data of Eq. ([1]). We then calculate the cross section for 
the DM-nucleus elastic scatterings and compare it with the current experimental bounds. 

3.1 Relic abundance 



The DM particle x interacts with a SM charged lepton through the coupling yi. When the 
temperature of the radiation-dominated Universe cools down below T ~ m^, x decouples from 
the thermal bath and its abundance freezes out. The decoupling process of x is regulated by 
the Boltzmann equation [8]. Let s, n^ and rf^, respectively, denote the entropy density, the 
number density of x ^i-iid the number density if x were to stay in thermal equilibrium. We 
further define the following quantities: 
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where g^, denotes the effective thermal degrees of freedom during the decoupling process of Xi 
and Mp = 2.4 x 10^^ GeV is the reduced Planck mass. For m^ < 300 GeV, we have g^, = 86.25. 
For z >3, n'^{z) can be written as 



n':''(z) 



9x 



^^2 ^ 



(10) 



where g^ = 2 counts the physical degree of freedom of the Majorana fermion x- 

We consider the DM pair annihilation process X X ^ ^j Cj, for which the thermally averaged 
cross section times velocity is expressed as [9] 
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^^') = 9i , eq. ... TTTT^ / ds. --m^i^W^ |-M(^)r , (H 



where |A^(s)p denotes the scattering rate per co-moving volume for the process and is given 
by 
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(12) 



The Boltzmann equation takes the form: 

dY{z) z {av) s 



{Y\z)-Yl{z)). (13) 



We numerically solve the Boltzmann equation for various values of y and ms- In our analysis, 
the mass of x is taken so as to satisfy the relations, m^ = ms — 50 GeV and m,^ = m^s — 20 
GeV, for simplicity. We evaluate the relic abundance of x as Y{z = oo). In Fig. 1, we show 
contour plots of the relic abundance of x oii the (m^, y^)-plane. Regions between two (red) 
lines correspond to the parameter regions where the resultant DM relic abundance fits with the 
observed values in Eq. ([1]). From Fig. 1, we can see that the mass difference between x £^nd S 
does not significantly alter the relic abundance as long as the difference is within 0(10 GeV). 
For the left panel in Fig. 1, we have found that the following relation must be satisfied in order 
to reproduce the observed DM abundance: 

2/2 ~ a{ms/GeV - 100) + b , (14) 

with 

0.0038 <a< 0.0040, 0.34 <b< 0.37 . (15) 
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Figure 1: Contour plot for the relic abundance of the DM candidate x o^i the (m^, |/^)-plane. 
We have taken m^ = ms — 50 GeV (m^ = ms — 20 GeV) in the left (right) panel. The region 
between two (red) lines corresponds to the range of the observed DM abundance in Eq. ([T]). 

3.2 Direct Dark Matter Detection 

The DM particle x couples to SM quarks at the one- loop level through the Feynman diagrams 
in Fig. 3u. We calculate effective couplings between x ci-nd a SM quark. The Higgs boson 
exchange leads to 
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where Q denotes a SM quark and x = mj'^/m?^. An effective coupling that arises from the 
Z-boson exchange is given by 
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where m^i denotes the mass of the SM lepton, gzi represents the coupling between the Z-boson 
and the SM lepton, and gzQ/^ denotes the axial part of the coupling between Z-boson and the 
quark Q. In deriving Eqs. f fT6|) and flT7|) . we have taken the limit m^i <C m^. 

Since the coupling between x ^i-iid a SM quark induced by the Z-boson exchange is of 
pseudo-scalar type, its contribution to DM-nucleus elastic scatterings is spin- dependent and 



^ Since x is a Majorana ferniion, photon exchange does not contribute, as it only gives rise to a magnetic 
dipole moment couphng. 
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Figure 2: Feynman diagrams contributing to the interaction between x ^^id a SM quark. The 
left diagram represents the Higgs boson exchange and the right one the Z-boson exchange. 

is proportional to a DM velocity squared in the Earth frame. Therefore this contribution is 
negligibly small. On the other hand, the coupling induced by the Higgs boson exchange diagram 
is of scalar-type and hence spin-independent. This contribution to the cross section of DM- 
nucleus elastic scatterings can be accessible by direct DM detection experiments if a value of 
K is large enough. In the following, we only consider the contribution from the Higgs boson 
exchange diagram. 

Let us calculate the cross section for DM-nucleon elastic scattering. We first estimate the 
following matrix elements: 

(iV| mggg |iV) = m^/^Q (g = u,ci,s,c,6,t) , (18) 

where iV represents a proton, p, or a neutron, n. A recent lattice calculation [10] reports that 

/^. + /^. = /t. + /™ =^ 0.056, |/^J<0.08. (19) 

We hereafter adopt these values, but we set /|^ = to make our analysis conservative. Opera- 
tors involving heavy quarks, c, 6, t, give rise to an effective coupling to gluons through a triangle 
diagram [11], and in heavy quark limit we have 



(iV| m^QQ \N) ~ (iV| (- Y^«.) G^uG^" \^) ^r Q = c, 6, t. (20) 

Using the trace of the QCD energy momentum tensor given by 

^M= E ^qQQ-^-^g^^g^^^ Y. ^qQQ-'^-^g.^g^^ , (21) 

Q=u,d,s,c,b,t Q=u,d,s 

we obtain 

(ATI J2 ^qQQ\N) + {N\ (^) G,.G^'''\N) . (22) 
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Figure 3: The spin-independent cross section for DM-nucIeon elastic scattering as a function 
of the DM mass, m^^,. Here we have fixed k = 3. The charged scalar mass, ms-, is related to the 
DM mass as m^ = ms — 50 GeV (lower two sohd lines) and m^ = ms — 20 GeV (upper two 
solid lines), respectively. The coupling constant y is taken to satisfy the relation of Eq. (!32|) so 
that the model reproduces the observed DM abundance. For each pair of degenerate lines, the 
upper one corresponds to (a = 0.0040, b = 0.37), while the lower one to (a = 0.0038, b = 0.34) 
in eq. dHl). The current experimental bound reported by the XENON 100 experiment ^ is 
shown by the dashed hne. 



We thus find for Q = c,b,t, 
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Using above expressions, the spin-independent cross section for the elastic scattering of x with 
a nucleon is given by 
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(24) 



where rriN ~ 1 GeV is the mass of a nucleon (proton or neutron), and x = m^/m^g. 

We evaluate o"si for various values of m^ = 7725 — 50 GeV and m^ = 7715—20 GeV, respectively. 
Here, y is given as a function of ms by using Eq. (TT4|) so as to reproduce the observed dark 
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matter relic abundance, and we fix the Higgs boson mass as nih = 125 GeV. In Fig. 3, we show 
(Tsi as a function of rriy. for k = 3, along with the current upper bound reported by the XENON 
100 experiment [3]. The resultant cross section is found to be (far) below the upper bound. 
Note that o"si scales as k^. If we allow a value of \k\ as large as its naive upper bound, \k\ < Air, 
from the view point of perturbation theory, the spin-independent cross section can be close 
to the current limit. The sensitivity of the spin-independent cross section in the future direct 
DM detection experiments, such as the XENON IT [12], is expected to reach (2 — 8) x 10~^^ 
cm~^, and hence our scenario can be tested. As we will discuss in the next sections, the spin- 
independent cross section has a correlation with a deviation of the Higgs-to-diphoton branching 
ratio from the SM predicted value. We will see that the deviation becomes larger as \k,\ is raised, 
equivalently, the spin-independent cross section becomes larger. 

4 Anomalous Higgs coupling to diphoton 

Although the Higgs boson discovered at the LHC experiments [5] has shown its properties 
mostly consistent with those expected in the SM, there exist some deviations from the SM 
expectations, namely, in its signal strength of the diphoton decay mode. Since the Higgs-to- 
diphoton coupling arises at the quantum level even in the SM, there is a good chance for certain 
new physics effects to alter the coupling. This has motivated many recent studies [I3l [H]. For 
studies in this direction before the Higgs boson discovery, see, for example, |15] . 
In our model, the charged scalar 5* couples to the Higgs boson through the term: 

C D -k{S^S){H^H). (25) 

Hence a one-loop diagram involving S gives a new contribution to the effective Higgs-to- 
diphoton coupling. The amplitude of the loop diagram of 5* is given by [16] 

As = C ^^\ [-XS + f{xs)] , (26) 

rrig Xg 

while the amplitudes in the SM via the top quark loop and the VT-boson loop are, respectively, 
given by 

At = C N, Q)\ {xt + [xt - l)/(xi)] , (27) 

Aw = -C ^[2xl, + 3xw + 3(2xw-l)fixw)] , (28) 



^w 



where C is a common constant, v = 174 GeV, A''^ = 3, Qf = 2/3, xs = m\/4m'^g, Xt = m1/4mf, 
and xw = f^'h/^''^w- "^^^ function /(x) is defined as 

f{x) = arcsin2(v^) (29) 



for a: > 1. 

We define the ratio of the branching fraction Br(/i — )■ 77) over its SM value as 
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(30) 



Note that as long as \As\ < \At + Aw\, As interferes constructively (destructively) with the 
SM contributions when k > (k < 0). If \k,\ is large and/or the scalar S is light, l^^l can 
cause a sizable deviation of the ratio from the SM prediction r^^ = 1. 

5 Implications to future experiments 

In our analysis, only three parameters {y, k, and ms) are involved with the fixed mass difference 
of m^ = ms — 50 GeV or m^ = ms — 20 GeV. In order to reproduce the measured relic 
abundance, we have found the relation between y and m,s in Eq. (TT4|) . by which only two 
parameters (k and m,s) are left free. As a result, both the spin-independent cross section in 
Eq. flM|) and the ratio of the Higgs-to-diphoton branching ratio in Eq. fl30l) can be obtained as 
a function of the two free parameters. Therefore, for a fixed m,s value, there is a correlation 
between o"is and r^^ through k. 

Fig. 4 shows the correlation for various values of m,s with the fixed mass difference of 
m^ = ms — 50 GeV. As \n\ is raised, the resultant spin- independent cross section increases 
and the deviation of r^^ from the SM value r^^ = 1 becomes larger. In the Figure, we have 
imposed the condition \n\ < 47r as a naive perturbative bound, which gives the upper bound on 
the cross section and the deviation of r^^ for a fixed m^ value. We have found that in order for 
the spin-independent cross section to be within the reach of the future experimental sensitivity 
(Tis ^ 10""^^ cm~^, the deviation of r^^ from the SM prediction should remain sizable. Therefore, 
the precise measurement of Higgs boson decay to diphoton at the LHC and the future direct 
DM detection experiments are complementary to reveal our scenario. 

6 Case for Dirac fermion dark matter 

In general, we can also consider a model similar to the one described in Sec. 2 except that x is 
a Dirac fermion. In this case, we introduce a global U(l) symmetry, instead of the Z2 parity in 
Table 1. We may assign U(l) charge -|-1 (—1) for the Dirac fermion DM (the complex scalar 
x). The basic Lagrangian remains the same as Eq. ([5]) (except for the correct normalization 
factor for the Dirac fermion). In this section, we will investigate phenomenology for the Dirac 
fermion case. 
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Figure 4: Correlation between r^^ and o"si through |fi:| < An with various ms values. The thick 
solid line, thin solid line, dashed line, dot-dashed line and dotted line, respectively, correspond 
to the charged scalar mass of ms = 100 GeV, 150 GeV, 200 GeV, 300 GeV and 400 GeV. The 
DM mass, m,^, is taken as m^ = ms — 50 GeV. 

Let us first evaluate the relic abundance of x i^i the same way as in Subsection 2.1 but for 
the Dirac fermion, |A^(s)p is replaced by 
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(31) 



with g^^ = 4. Numerically solving the Boltzmann equation, we find the DM relic abundance as a 
function of m,s and y"^. Corresponding to Fig. 1, we show our result in Fig. 5, for m^ = m^s — 50 
GeV. As in the case of Majorana fermion DM, the difference between ?7i^ and ms does not 
significantly alter our results as long as the difference is within 0(10) GeV. From Fig. 5, we 
find that the relation 
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(32) 
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Figure 5: Contour plot for the relic abundance of the DM candidate x on the {ms, |/^)-plane. 
Here we have taken m^ = rris — 50 GeV. The region between two (red) lines corresponds to the 
range of the observed DM abundance in Eq. ([T]). 



with 



0.0025 < a < 0.0027, 0.29 <b< 0.31 



(33) 



in order to reproduce the range of the observed relic abundance in Eq. ([T]). 

As in the Majorana case, the DM couples to SM quarks at the one-loop level. In addition 
to the diagrams in Fig. 2, there is a diagram mediated by photon, shown in Fig. 6, through the 
magnetic dipole moment induced at one- loop level. In fact, this new diagram gives a dominant 
contribution to the spin- independent cross section for the DM elastic scattering with nucleuC- 
The effective coupling induced by the magnetic dipole moment is found to be 
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where g'^ denotes the momentum transfer from x to a quark (Q), g^i {g-yo) denotes the QED 
coupling of the SM lepton (the SM quark), and x = m^/m^. 



For general discussions for a DM particle with electric/magnetic dipole moments, see Ref. |17) . 
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Figure 6: Feynman diagram that gives a dominant contribution to the interaction between x 
and a SM quark when x is a Dirac fermion. 

Since the vector interactions with each quarks coherently contribute to the total cross sec- 
tion, we obtain the differential cross section for x-nucleus elastic scattering (in the center-of- 
mass frame) as 

^2 - 2mr 
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(35) 



where rriN and Z respectively denote the mass and the atomic number of a nucleus. The 
divergence at 6* = is cut off by the low energy threshold for a recoiling nucleus. To make a 
conservative comparison of the cross section in our model with the bounds reported by direct 
DM detection experimentqj, we replace cos 6* — )■ —1 in the right hand side of Eq. (l35ll . Then 
the conservatively-estimated total cross section for ^-nucleus elastic scattering is given by 
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2 H r 

X 



IGvr' 



-,y 



m^ X + log(l 



x] 



ml. 



2x2 



.(36) 



(m7V5 + m^)2 

We finally translate the DM-nucleus elastic scattering cross section into the spin-independent 
cross section for DM-nucleon elastic scattering, {(Tsi) as follows: 



0"SI 



(^X-NS X 



mil {f^NS + m^f ]_ 

^lis {mN + m^Y ^2 ' 



(37) 



^ Direct DM detection experiments actually put a bound on the normalization for the distribution of the 
event rate per recoil energy, dR/dEr {R denotes the event rate and Er the recoil energy), with the assumption 
that dR/dEr oc cxp{~ AEr), where A is some known constant. In our case, however, we have dR/dEr ex 
(1 + B / Er) exp{— AEr) with a constant B because of the term 1/(1 — cos9) in Eq. (|35|) . Thus, we cannot 
naively compare the elastic scattering cross section in our model with an experimental bound. All we can do 
is to make a conservative estimate on the cross section by substituting —1 into cos 9, and compare it with the 
experimental bound. 
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Figure 7: The conservatively-estimated spin-independent cross section for DM-nucleon elastic 
scattering (o"si) as a function of the DM mass (rn^). The charged scalar mass, m^, is related to 
the DM mass as m^ = rris — 50 GeV (lower two solid lines) and m^ = ms — 20 GeV (upper two 
solid lines). The coupling constant y is taken to satisfy the relation in Eq. ( 132|) so that the model 
reproduces the observed DM abundance. For each pair of degenerate solid lines, the upper line 
corresponds to (a = 0.0027, b = 0.31) in Eq. §^, while the lower one to (a = 0.0025, b = 0.29) 
For comparison, the current experimental bound reported by the XENON 100 experiment [3] 
is shown by the dashed line. 

where A denotes the mass number of the nucleus. For the XENON 100 experiment, we take 
Z = 54 and A = 131.3. 

In Fig. 7, we show the conservatively-estimated spin-independent cross section for x-nucleon 
elastic scattering, o"gi, as a function of m^ which is taken to be m^ = ms — 50 GeV and 
m^ = ms — 20 GeV, respectively. Here, y is fixed to satisfy the relation of Eq. ( 132|) . We 
see that the resultant cross sections are about an order of magnitude smaller than the current 
upper bound and will be tested in future direct DM detection experiments. 
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